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INTRODUCTION 


The  Ml  19  A1  howitzer  is  a  lightweight  howitzer  used  by  the  U.S.  Army.  Since  its 
introduction  to  the  U.S.  Army  in  1987,  the  Ml  19  howitzer  has  undergone  numerous  design 
changes  to  improve  maintainability,  functionality,  and  compatibility.  It  is  currently  being  used  in 
the  Middle  East  by  regular  Army  and  by  the  National  Guard  (fig.  1). 


Figure  1 

Soldiers  train  on  firing  points  with  an  Ml  19  howitzer  (ref.  1) 

This  report  describes  modeling  and  simulation  of  the  thermal  effects  of  a  new  power 
enclosure  for  the  Ml  19  howitzer.  The  power  enclosure  will  increase  functionality  of  the  mobile 
unit.  The  enclosure  consists  of  a  sealed  box  containing  batteries,  power  control  circuitry,  and 
other  assorted  electronics.  The  enclosure  regulates  and  supplies  power  to  electronics.  To 
maintain  reliability  of  the  power  supply  unit,  each  component  must  remain  below  its  design 
temperature  during  expected  operating  conditions.  For  this  initial  study,  it  was  assumed  that  the 
gun  vibration  and  temperature  rise  in  the  gun  tube  would  not  affect  the  power  enclosure  (refs.  2 
through  4). 

The  thermal  analysis  was  done  on  the  enclosure  using  the  extreme  operating  conditions. 
The  operating  conditions  account  for  extreme  solar  radiation,  elevated  ambient  temperature, 
and  heat  generation  of  components.  The  system  was  modeled  using  Computational  Fluid 
Dynamics  (CFD).  This  computer  model  was  used  to  solve  the  heat  transfer  problem  including 
conduction,  convection,  and  radiation.  This  approach  allows  for  changes  to  the  parameters  of 
the  model  as  needed  for  redesign  or  upgrades. 

Two  thermal  load  cases  were  evaluated.  First,  the  enclosure  and  the  parts  inside  the 
enclosure  were  modeled  using  conduction  with  a  prescribed  convection  on  the  outside  of  the 
enclosure.  The  second  analysis  was  built  upon  the  first  analysis.  Free  convection  on  the  outer 
surfaces  of  the  enclosure  and  solar  radiation  were  added  to  the  thermal  loads.  The  goal  of  both 
of  these  analyses  was  to  find  the  temperature  of  various  components  once  the  system  reached 
a  steady  state  thermal  condition  during  operation. 
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METHODS  ASSUMPTIONS  AND  PROCEDURES 


Geometry 

For  this  analysis,  geometry  was  imported  from  Pro/E  (ref.  5)  part  files  for  the  enclosure 
and  the  internal  components.  The  model  as  provided  had  details  that  would  not  affect  the 
results  and  make  the  model  overly  complex.  To  simplify  the  model,  only  major  components 
were  considered.  The  analysis  included  the  box,  battery,  main  relay,  heat  sink,  circuit  card 
board,  charger  board,  burst  board,  and  the  DC  filter  board.  The  inside  of  the  enclosure  was 
cluttered  with  wires  and  cables.  None  of  these  were  included  in  this  simplified  model.  ProE  files 
were  converted  to  CFD  files  using  the  software  package  Gambit  (ref  6).  Figure  2  shows  the 
Pro/E  assembly. 


Relay 


Com  muni  cation 


DC  Filter 
Circuit  Card 


Burst  Supply 
Circuit  Card 


Battery 


Charger 
Circuit  Card 


Figure  2 

Pro/E  geometry  of  power  supply  package 

Model  Setup  and  Assumptions 

The  geometry  files  were  imported  to  the  CFD  software  package,  STAR-CCM+  v 
3.06.006  (ref  7).  The  imported  model  included  the  enclosure  box,  battery,  main  relay,  heat  sink, 
circuit  card  board,  charger  board,  burst  board,  and  the  DC  filter  board. 

Once  the  geometry  was  in  the  software  package  STAR-CCM+,  it  was  further  refined. 
This  refinement  reduced  the  number  of  contact  interfaces  without  any  adverse  impact  on  the 
results  of  the  simulation.  To  model  the  transfer  of  heat  from  one  body  to  another,  contact 
interfaces  needed  to  be  defined  between  surface  faces  that  were  touching  each  other. 
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All  of  the  components  were  modeled  as  3-D  stationary  solids  using  the  coupled  solid 
energy  solver  and  having  constant  density.  This  included  the  air  within  the  enclosure.  It  was 
assumed  that  the  wire,  cables,  and  other  internal  blockages  would  allow  little  natural  convection 
within  the  enclosure.  The  air  was  modeled  as  a  stationary  solid  that  could  only  conduct  heat. 

Materials 

This  analysis  consisted  of  many  different  parts  each  with  their  own  thermal  properties. 
For  simplicity  and  due  to  lack  of  data,  everything  was  assumed  to  be  one  of  seven  possible 
materials.  For  each  material  the  density,  thermal  conductivity,  and  heat  capacity  were  needed. 
Some  material  parameters  were  found  in  the  open  literature  and  others  were  estimated  from 
available  references  (refs.  8  and  9). 

The  components  were  assumed  to  be  in  one  of  seven  material  groups.  Table  1  shows 
the  material  properties  for  air,  printed  circuit  board  (PCB),  aluminum,  battery,  electronics,  large 
power  chips  (power  converters),  and  small  power  chips  (power  filters). 


Table  1 

Material  assumptions 


Material 

Density  (ka/mA3) 

Thermal  conductivity 
(W/mK) 

Heat  capacity 
(W/m-K) 

Air 

1.18415 

0.0257 

0.026 

PCB 

1900 

23 

1200 

Aluminum 

2800 

193 

880 

Battery 

534 

84.7 

3582 

Electronics 

3600 

17.17 

795.5 

Large  power  chips 

950 

165 

120 

Small  power  chips 

950 

100 

83 

Interfaces 

The  areas  where  two  parts  touch  are  defined  as  contact  interfaces;  STAR  requires  a 
thermal  contact  resistance  between  two  contact  surfaces.  Contact  resistance  is  a  function  of 
surface  finish,  contact  pressure,  and  interface  material  (ref.  8).  Since  this  data  was  not  available 
for  the  components,  an  equivalent  thermal  contact  resistance  was  found  using  an  air  gap  size 
associated  with  standard  surface  finishes.  For  the  interfaces  where  thermal  compound  was 
applied,  that  compound’s  thermal  conductivity  was  used  in  place  of  air’s.  These  values  were 
compared  with  the  available  data  in  literature  and  found  to  be  within  the  correct  ranges  (ref.  8). 
Interface  assumptions  are  listed  in  table  2. 
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Table  2 

Calculated  contact  interface  resistance 


Thermal  contact  resistance 


Interface 


(mA2-K/W) 


AL-PCB 

AL-AL  (thermal  grease) 
AL-power  chips 

AL-power  chips  (thermal  grease) 
AL-electronics 
PCB-power  chips 
PCB-electronics 


0.000389 

0.000083 

0.003891 

0.000033 

0.0389 


0.019 

0.019 


Operating  Conditions 

The  system  is  designed  to  operate  reliably  in  an  area  with  an  ambient  temperature  of 
130°F  (327.6  K).  In  addition  to  the  ambient  requirement,  the  power  supply  unit  must  operate 
with  a  solar  load  of  1 1 20  W/mA2.  Figure  3  shows  the  view  of  the  box  that  the  source  of  the  solar 
radiation  sees.  The  inside  of  the  box  is  cluttered  with  cables,  wiring,  and  other  small  parts  that 
were  not  modeled.  This  clutter  may  block  airflow  at  critical  locations  within  the  box.  To  take  this 
blockage  and  lack  of  internal  air  flow  into  account,  the  air  inside  the  box  was  modeled  such  that 
it  only  factored  into  the  heat  transfer  calculations  as  allowing  heat  transfer  through  conduction 
only. 


Figure  3 

Solar  radiation  view  of  the  enclosure 


The  heat  generation  for  the  various  components  was  modeled  as  a  constant  value  in 
time.  These  values  were  calculated  from  the  known  or  estimated  heat  production  rates  and  are 
expressed  in  W/mA3.  Table  3  shows  heat  production  and  heat  produced  per  unit  volume  for 
each  component. 
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Table  3 

Heat  generation  by  internal  components  in  the  power  supply  unit 


Component 

QiW} 

W/mA3 

Burst  chip 

23 

277108 

DC  filter  big  chips 

7 

176322 

DC  filter  small  chips 

8 

390224 

Relay 

0.5 

714 

Battery 

1 

131 

Circuit  card 

2 

6779 

Charger  main  chip 

33 

397590 

Charger  chip  A 

0.85 

798122 

Charger  chip  B 

26 

1150442 

Charger  chip  C 

2.6 

1135371 

Charger  chip  D 

0.4 

197044 

Charger  chip  E 

2.5 

1288660 

Conduction  only,  Analysis  1 

The  first  simulation  was  done  without  solar  radiation  or  natural  convection.  However, 
since  natural  convection  would  be  an  important  method  of  heat  transfer  out  of  the  system,  it  had 
to  be  taken  into  account.  This  was  accomplished  by  specifying  an  ambient  temperature  (130°F) 
and  a  conservative  convective  heat  transfer  coefficient  (5  W/mA2)  to  all  of  the  outer  surfaces. 
These  outer  surfaces  would  normally  be  in  ‘contact’  with  the  atmosphere.  This  provided  a 
relatively  quick  analysis  that  formed  the  ground  work  and  a  bench  mark  for  further  analyses. 

Solar  Radiation  and  Natural  Convection,  Analysis  2 

The  model  was  improved  by  including  solar  radiation  and  natural  convection.  The 
convection  analysis  included  some  updated  heat  production  rates  (table  4). 

Table  4 

Updated  heat  production  rates  for  selected  components,  analysis  2 


Component 

Q(W) 

W/mA3 

DC  filter  small  chips 

9.25 

4512195 

Charger  chip  D 

0.63 

310345 

Charger  chip  E 

3.85 

1984536 

The  power  enclosure  was  designed  to  withstand  relatively  high  solar  loads  (ref.  10).  The 
direct  solar  load  was  modeled  as  1120  W/mA2  and  the  diffuse  solar  load  as  840  W/mA2  (75%  of 
the  direct  solar  load).  These  values  assumed  that  the  sun  was  in  the  “worst”  position  in  the  sky. 
The  direct  solar  load  was  applied  to  the  sides  of  the  box  with  the  components  of  interest.  The 
direct  solar  load  was  also  applied  to  the  top  of  the  box.  The  diffuse  solar  load  was  applied  to  the 
rest  of  the  exterior  surfaces,  modeling  the  scatting  effect  of  the  atmosphere  and  the  reflections 
from  objects  that  could  be  nearby. 
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Since  natural  convection  is  a  buoyancy-driven  flow,  the  gravity  model  in  STAR  also  had 
to  be  activated.  The  outer  surfaces  of  the  box  and  heat  sink  had  their  boundary  conditions 
modified  to  be  compatible  with  the  natural  convection.  To  be  able  to  model  the  natural 
convection,  it  was  necessary  to  add  a  volume  of  air  around  the  box.  This  did  not  change  the 
internal  components. 


RESULTS 


Conduction  Only 

Figure  4  shows  the  results  of  the  conduction-only  analysis.  With  conduction  only,  Diode 
B  and  Relay  2  exceeded  failure  temperatures.  Other  component  temperatures  neared  the  failure 
condition,  but  did  not  exceed  design  values. 


Temperature  (C) 
131.00  15400 


177.00 


Figure  4 

Temperature  results,  conduction  only 
Conduction,  Solar  Radiation,  and  Natural  Convection 

Figure  5  shows  the  result  with  radiation  and  convection  included.  Figure  6  shows  the 
components  that  exceed  their  design  temperature  values.  With  solar  radiation  and  natural 
convection,  more  than  half  of  the  components  would  be  expected  to  exceed  their  design 
temperatures  at  extreme  conditions. 
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Figure  5 

Temperature  distribution  with  convection,  conduction  and  radiation  included 


Converters- 


Resistors: 


Failure  at  1 00C  Failure  at  275C 


Filters. 

Failure  at  100C 


- raiiureai^oo  Diode  As 


Diode  Bs: 
Failure  at  HOC 


Relays- 

Failure  at  140C 


Figure  6 

Summary  of  failed  components  based  on  thermal  analysis 
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Table  5  compares  the  results  with  and  without  solar  radiation  and  convection.  Solar 
radiation  and  convection  increased  the  temperature  in  key  components  by  about  20%. 

Table  5 

Component  temperatures 


Failure  temperature 

Conduction 

w/Solar  radiation  and  natural 

f°C) 

only  (°C) 

convection  f°C) 

Converter  -  burst 

100 

98.4 

112.2 

Converter -charger 

100 

94.1 

105.5 

Filter 

100 

97.9 

117.1 

Diode  A 

160 

149.1 

159.7 

Diode  A  2 

160 

147.7 

158.1 

Diode  B 

110 

113.1 

140.7 

Resistor 

275 

129.3 

140.1 

Relay  1 

140 

103.1 

119.1 

Relay  2 

140 

153.9 

198.1 

CONCLUSIONS 

The  current  design  is  unable  to  keep  the  components  within  the  allowable  operating 
temperature.  As  expected,  the  added  solar  load  in  the  second  simulation  resulted  in  higher 
temperatures  and  additional  failed  components. 

Several  design  improvements  could  be  suggested.  The  system  does  have  an  external 
heat  sink  to  help  dissipate  the  heat  that  was  included  in  this  model.  This  heat  sink  could  be 
optimized  to  provide  greater  heat  dissipation.  In  addition,  it  may  be  possible  to  install  an 
additional  heat  sink  to  dissipate  heat  from  the  side  with  the  filter  and  burst  chips.  A  third 
suggestion  is  to  replace  components  with  similar  components  that  generate  less  heat  and/or 
function  at  higher  temperatures. 


8 


REFERENCES 


1.  Olive-Drab.com;  ttp://www. olive-drab.com/idphoto/id_photos_m1 19_105mm.php;  22 
May  2008. 

2.  Carlucci,  P.  and  Potucek,  J.;  “Transient  Thermal  Response  of  a  Projectile  Misfired  in  a 
Mortar  Tube;"  Technical  Report  ARMET-TR-08017;  U.S.  Army  Armament  Research, 
Development  and  Engineering  Center,  Picatinny  Arsenal,  NJ;  September  2008 

3.  Cordes,  J.  A.  and  Reinhardt,  L.  Reinhardt,  P.E.;  The  Rise  and  Fall  of  Temperature  in  a 
Gun  Tube  Due  to  Multiple  Gun-Firings;  21st  International  Symposium  on  Ballistics, 
Adelaide,  South  Australia,  pp.  651-654,  19-23  April  2004. 

4.  Carlucci,  D.E.  and  Jacobson,  S.S.;  Ballistics  Theory  and  Design  of  Guns  and 
Ammunition.  CRC  Press,  NYC,  NY,  2008. 

5.  Pro-E 

6.  Gambit 

7.  Star-ccm 

8.  Mills,  A.F.;  Heat  Transfer;  2nd  Ed.,  Prentice-Hall,  Upper  Saddle  River,  NJ;  1999. 

9.  Pecht,  Michael;  Handbook  of  Electronic  Package  Design;  Marcel  Dekker,  Inc.  New  York, 
NY;  1991. 

10.  Army  Regulation  70-38;  Research,  Development,  Test  and  Evaluation  of  Materiel  for 
Extreme  Climatic  Conditions;  September  1979. 


BIBLIOGRAPHY 

Ham,  Jeff;  Email  communication  concerning  thermal  properties  of  electronic  components',  Vicor 
25  Frontage  Rd,  Andover  MA;  February  2009. 

LeFante,  Michael;  Email  communication  concerning  geometry,  thermal  properties,  and  setup 
conditions  for  the  model;  U.S.  Army  Armament  Research,  Development  and 
Engineering  Center,  Picatinny  Arsenal,  NJ;  February  2009. 


9 


DISTRIBUTION  LIST 


U.S.  ArmyARDEC 
ATTN:  RDAR-EMK 
RDAR-GC 
RDAR-MEF-E  (10) 

Plcatinny  Arsenal,  NJ,  07806-5000 

Defense  Technical  Information  Center  (DTIC) 

ATTN:  Accessions  Division 
8725  John  J.  Kingman  Road,  Ste  0944 
Fort  Belvoir,  VA  22060-6218 

Commander 

Solder  and  Biological/Chemical  Command 
ATTN:  AMSSB-CII,  Library 

Aberdeen  Proving  Ground,  MD  21005-5066 

Chief 

Benet  Weapons  Laboratory,  WSEC 

U.S.  Army  Research,  Development  and  Engineering  Command 
Armament  Research,  Development  and  Engineering  Center 
ATTN:  RDAR-WSB 
Watervliet,  NY  12189-5000 

Director 

U.S.  Army  TRADOC  Analysis  Center-WSMR 

ATTN:  ATRC-WSS-R 

White  Sands  Missile  Range,  NM  88002 

Chemical  Propulsion  Information  Agency 
ATTN:  Accessions 
10630  Little  Patuxent  Parkway,  Suite  202 
Columbia,  MD  21044-3204 

GIDEP  Operations  Center 
P.O.  Box  8000 
Corona,  VA  91718-8000 

Commander 

Air  Force  Research  Laboratory 
ATTN:  Technical  Library 
2977  P  St.,  STE  6 

Wright  Patterson  AFB,  OH  45433-7739 
Commander 

Air  Force  Research  Laboratory 
101  W.  Eglin  Blvd 
EglinAFB,  FL  32542-6810 
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Commander 

Naval  Surface  Warfare  Center 
Dahlgren  Division  Code  G06 
ATTN:  Tech  Library,  Code  323  (6) 
17320  Dahlgren  Rd. 

Dahlgren,  VA  22448 


